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Abstract We investigate the low-temperature behavior of the ratio between the 
shear viscosity 77 and the entropy density 5 in the unitary Fermi gas by using a 
model based on the zero-temperature spectra of both bosonic collective modes and 
fermonic single-particle excitations. Our theoretical curve of r]/s as a function of 
the temperature T is in qualitative agreement with the experimental data of trapped 
ultracold 6 Li atomic gases. We find the minimum value r}/s ~ 0.44 (in units of 
h/kg) at the temperature T/Tp ^ 0.27, with Tp the Fermi temperature. 
CT 1 ! 
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1 Introduction 

Strongly interacting quantum many-body systems like Helium 4, the quark-gluon 
plasma and the unitary Fermi gas share a common feature: an extremely low 
viscosity hydrodynamics. 12 These quite different many-body systems show a 
ratio of shear viscosity 77 to entropy density 5 which is not too far from the 
lower bound r\/s = h/(4nkB) predicted for a "perfect fluid" by using the anti- 
deSitter/conformal field theory (AdS/CFT) duality between certain strongly cou- 

lO ■ pled field theories in d = 4 space-time dimensions and weakly coupled string the- 

ory in d = 10. 3 As discussed in a recent review, 7 theoretical predictions 5 ' 6 of the 
viscosity-entropy ratio for dilute and ultracold Fermi atoms in the unitary regime, 

(^. ■ where the s-wave inter-atomic scattering length cif diverges, are not in good agree- 

ment with the experimental data of the viscosity-entropy ratio measured in the 6 Li 
atomic gas.- 

In this paper we study the low-temperature behavior of r\ /s by using a recent 
heuristic analysis of the shear viscosity 8 and a thermodynamical model 9 of the 
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unitary Fermi gas based on zero-temperature elementary excitations. We show that 
our theoretical curve for 77 /s as a function of the temperature T is in qualitative 
good agreement with the experimental data of trapped ultracold 6 Li atomic gases. 
In particular, we find the minimum value r\/s ~ 0.44 (in units of h/kg) at the 
temperature T /Tp a 0.27, with Tp the Fermi temperature. Both the value and the 
position of this minimum are fully compatible with the most recent experimental 
determinations .2A 

In the first part of this paper we briefly review our thermodynamical model 9 
of the unitary Fermi gas comparing it with experimental data 10 and Monte Carlo 
simulations 1 1 . In the second part we adopt the analysis of How and LeClair 8 for 
the shear viscosity and derive from it and from our thermodynamical model 9 the 
viscosity-entropy ratio 77 /s. We then compare our curve of r)/svsT with available 
experimental data 4 and proposed theories. 5 i 6 i 7 i 8 

2 Elementary excitations of the unitary Fermi gas 

For any many-body system the weakly excited states, the so-called elementary 
excitations, can be treated as excitations of an ideal gas. 12 i 13 In general, these el- 
ementary excitations are the result of collective interactions of the particles of the 
system, and therefore pertain to the system as a whole and not to its separate par- 
ticles. 1213 For the unitary Fermi gas the mean-field extended BCS theory predicts 
the existence of fermionic single-particle elementary excitations characterized by 
an energy gap A . 14 The inclusion of beyond-mean-field effects, namely quantum 
fluctuations of the order parameter, gives rise to bosonic collective excitations,— 
which are density waves reducing to the Bogoliubov-Goldstone-Anderson mode 
in the limit of small momenta.— 

Our effective quantum Hamiltonian 9 of the uniform unitary Fermi gas with 
two equally-populated spin components is then assumed to be: 

H = E + Y,£coi(q) b^b q + £ £%,(» c+ a c pa , (1) 

q <T=t4. p 

where Eq is the ground-state energy, bt and b p are the bosonic creation and de- 
struction operators of a collective excitation of linear momentum q with energy 
£coi(q), while c+ and c po are the fermionic creation and destruction operators of 
a single-particle excitation of linear momentum p and spin a, with energy e sp (p). 
It is now well-established 14 that the ground-state energy Eq of the uniform 
unitary Fermi gas made of N atoms in a volume V is given by 

E = l^Ne P , (2) 

with I ~ 0.4 16 and where Ef = h 2 (37Z 2 n) 2 / 3 /(2m) is the Fermi energy of a non- 
interacting fermi gas with density n = N/V. 

The exact dispersion relation of elementary (collective and single-particle) ex- 
citations is not fully known. 14 In Ref. 15 we have found the dispersion relation of 
collective elementary excitations as 



.e,„ / U/i-\M< ? 2 + i ^ 4 , (3) 



where 



c\ = \l - vp , (4) 



is the zero-temperature first sound velocity, with vp = (h/m)(37t 2 n) 1 ^ 3 the Fermi 
velocity of a noninteracting Fermi gas. Notice that the term with k takes into ac- 
count the increase of kinetic energy due to spatial variations of the density- 



For the purposes of the present paper, by fixing £, = 0.42, i.e. the value given by 
the Monte Carlo prediction for a uniform gas of Astrakharchik et al.,— we find 
that the best agreement with Monte Carlo data is obtained with k = 0.25. 

The collective modes describe correctly only the low-energy density oscil- 
lations of the system while at higher energies one expects the appearence of 
fermionic single-particle excitations starting from the threshold above which there 
is the breaking of Cooper pairs. 14 ' 11 ' 24 At zero temperature these single-particle 
elementary excitations can be written as 



where £ is a parameter which takes into account the interaction between fermions 
(£ ~ 0.9 according to recent Monte Carlo results 24 ) with Ep the Fermi energy 
of the ideal Fermi gas. Aq is the zero-temperature gap parameter with 2Aq the 
minimal energy to break a Cooper pair. 14 Notice that the gap energy Aq of the 
unitary Fermi gas at zero-temperature has been calculated with Monte Carlo sim- 
ulations 24 ' 25 and found to be Aq — yep, with j— 0.45. 



3 Thermodynamics of the unitary Fermi gas 

At very low temperature the thermodynamic properties of the superfiuid unitary 
Fermi gas can be obtained from the collective spectrum and considering it as an 
ideal Bose gas of elementary excitations 12 with the bosonic distribution 

*® = <w = ^wL-i < (6) 

where (A) = Tr[Ae~ H l kBT }/Tr[e~ H l kBT } is the thermal average of the operator 
A with T the absolute temperature and kg is the Boltzmann constant. 26 As T 
increases also the fermionic single-particle excitations become important. Thus 
there is also the effect of an ideal Fermi gas of single-particle excitations with the 
fermionic distribution 

f F (p) = (e^Cpa) = e£sp{p)/kBT + l - (7) 

which is spin independent. 

The Helmholtz free energy F of any thermodynamic system is given by 

F = -k B T]n3f, (8) 



where 



3t = Tr[e-"/ kBT ] , 



(9) 

is the partition function of the system. 26 Using Eq. (QJ the free energy of our 
unitary Fermi gas can be written as F = Fq+F co i+F s „, where Fq is the free energy 
of the ground-state, F co j is the free energy of the bosonic collective excitations 
and F sp the free energy of fermionic single-particle excitations. The Helmholtz 
free energy Fq of the uniform ground state coincides with the zero-temperature 
internal energy Eq and is given by 



r$Ne F , 



(10) 



where N is the number of atoms of the uniform system in a volume V. The free 
energy F co i of the collective excitations is instead given by 12 



?col 
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-e c0 l(q)/(k B T) 



while the free energy F sp due to the single-particle excitations is 



(11) 



F sp = -2 k B T J^ln 



-es P (p)/(k B T) 



(12) 



where the factor 2 is due to the two spin components. As previously discussed, 
the total Helmholtz free energy F of the low-temperature unitary Fermi gas can 
be then written as Fq + F co i + F s „, namely 



N£ F <P[ — 
\T F 



(13) 



where <P(x) is a function of the scaled temperature x = T/Tp, with 7> = Sp/ks, 
given by 



3 3 

<P(x) = -$ + -x 

5 2 jo 
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l+e sp 



e sp {u)/x 



u du . 



(14) 



Here the discrete summations have been replaced by integrals, moreover we set 
e co i{u) = y/u 2 {Xu 2 +4£/3) and £ sp {u) = V V - Q 2 + 1 '■ 

From the Helmholtz free energy F we can immediately calculate the chemical 
potential /I, through 



dN 



T.V 



obtaining 



\x = e F 



3 \T F 



3 7> \T F 



(16) 



where <P'(x) 



d<P(x) 
dx 



and one recovers jIq = £, £p in the limit of zero-temperature. 
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Fig. 1 (Color online). Scaled internal energy E/{Ne F ) as a function of the scaled temperature 
T/Tp. Filled circles: Monte Carlo simulations 11 . Open squares with error bars: experimental 
data of Horikoshi et al. 10 . Solid line: our model, i.e. Eq. l l20t and Eq. i THT l. with X = 0.25. Solid 
line: our model, i.e. Eq. ( l20t and Eq. d!4l . with X = 0. Other zero-temperature parameters of 
elementary excitations: t, = 0.42, £ = 0.9, and y = 0.45. 



The entropy 5 is related to the free energy F by the formula 



2V,V 



from which we get 



s=-m*>(£) 



In addition, the internal energy E, given by 

E = F + TS . 



can be written explicitly as 



E=N£ F 



*lL)-L#(L 

T F T F \T F 



(17) 



(18) 



(19) 



(20) 



It is interesting to compare our model with other theoretical approaches and 
also with the available experimental data. In Fig. Q] we report the data of internal 
energy E obtained by Bulgac, Drut and Magierski ' ' with their Monte Carlo sim- 
ulations (filled circles) of the atomic unitary gas. We insert also the very recent 
experimental data of Horikoshi et al.— for the unitary Fermi gas of 6 Li atoms 
but extracted from the gas under harmonic confinement (filled squares with error 
bars). In the figure we include the results of our model, that is given by Eqs. J20I) 
and (fl4l) with both A = 0.25 (solid line) and A = (dashed line). The figure shows 



that in our model the gradient term, proportional to A, plays a marginal role up to 
T /Tp ~ 0.25. Above T/Tp ~ 0.25, however, our results with A = 0.25 shows a 
better agreement with both Monte Carlo and experimental data than those with 
A = 0. We stress that the gradient term is essential to describe accurately the zero- 
temperature surface effects of a trapped system, in particular with a small number 
of atoms, where the Thomas-Fermi (i.e. A = 0) approximation fails 15 . The value 
A = 0.25 gives the best fit of the Monte Carlo energy as a function of the particle 
number for £ = 0.42 (see Ref.-i^ for details). 

Our model does not show a phase transition. Nevertheless, the results of Fig.Q] 
strongly shows that our model works quite well not only in the superfluid regime, 
but also slightly above the critical temperature (T c /Tp ~ 0.15) suggested by two 
theoretical groups. 24 - 27 This finding is not fully surprising. In presence of a pseudo- 
gap region, the temperature-dependent gap A (T) of single-particle elementary ex- 
citations can be written as A(T) = A SC (T) +A pg (T), where As c (T) is the super- 
conducting gap and A pg (T) is the pseudogap. 14 At T c the superconducting gap 
A SC (T) goes to zero, i.e. A SC {T C ) = 0, but the pseudo-gap A pg (T) remains finite 
and it becomes zero only at the higher temperature T* . 14 For further details on the 
comparision between our model and other theories see Ref . 9 . 



4 Shear viscosity from thermodynamics 

A first principle calculation of the shear viscosity is beyond the scope of the 
present work and we adopt the heuristic analysis of How and LeClair 8 to write 
it in terms of the scaled free energy <P(x) and its first derivative <P'(x). The shear 
viscosity T] can be estimated by using the formula— 

ri = -nmvl m , (21) 

where n is the total number density of the fluid, m is the mass of each particle in 
the fluid, v is the average velocity of particles, and l m is the length of the mean 
free path. The mean free path is written as 

lm = — , (22) 

no 

where a is a suitable transport cross-section. 28 

For the unitary gas with two-spin-component fermions, the cross-section is 
given by 

An 

|ki-k 2 | 2 

where ki — k 2 is the relative wave number of two colliding fermions with opposite 
spin. 12 The average velocity v of fermions can be related to the relative wave 
number |kj — k 2 | by the formula 8 

v = ^2-<|k 1 -k 2 | 2 ) 1 /2. (24 ) 

m 
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Fig. 2 (Color online). Viscosity-entropy ratio r\/s (in units oihlkg) as a function of the scaled 
temperature T/Tp. Dot-dot-dashed line: the bound r\/s = \/{An) predicted by field theory mod- 
els inspired by string theory. 3 Filled circles with error bars: experimental data. 4,7 Dot-dashed 
line: low-temperature prediction of Ref. 5 . Dotted line: high-temperatue prediction of Ref. 6 . 
Solid line: our model with X = 0.25. Dashed line: our model with X = 0. Other zero-temperature 
parameters of elementary excitations: £, = 0.42, C, = 0.9, and y= 0.45. 



In fact, (|ki - k 2 | 2 ) = {k\ + k 2 - 2ki • k 2 ) = 2k 2 , because (ki • k 2 ) = and k 2 = 
(k 2 ) = (k 2 ), and finally v = (h/m)k. 2<} In this way the shear viscosity becomes 



n = 



3-3 

in v 
6nh 2 



(25) 



The average velocity v can be estimated by imposing that the average kinetic en- 
ergy is equal to the internal energy per particle, 8 namely 



1 2 E 

—mv = — 

2 N 



By using Eq. d25b with v given by Eq. 
viscosity reads 



71 = nh- (<P{ — ) - —#'( — ) 



(26) 
and E given by Eq. (120b . the shear 

3/2 



(27) 



Notice that for T — > 0, the viscosity 77 goes to a constant value because also a 
goes to a constant ~ kj- . This is in agreement with recent experimental results 
on the universal spin diffusion in a strongly interacting Fermi gas. 30 Finally, by 
considering Eq. ( 118b for the entropy of the unitary Fermi gas, we find that the 
viscosity-entropy ratio is given by 



1 

s 



~k B 2 



3/2 



*U: 



(28) 



where s = S/V is the entropy density, i.e. the entropy per unit of volume. This for- 
mula gives the viscosity-entropy ratio in terms of the scaled free energy <P(x) and 
its first derivative <P'(x). For T — > Eq. ( 128b gives r}/s = +°°. This divergence of 
T)/s is a consequence of Eqs. ( 125b and ( 126b which impose, as previously stressed, 
a small but finite viscosity 77 while the entropy density 5 goes to zero. 

In Fig. |2] we plot experimental data of the ratio J]/s (filled circles with error 
bars). These data have been obtained by the group of Thomas 4 from the damping 
of radial breathing mode of the atomic cloud, and then elaborated by Schafer and 
Chafin 7 with an energy-to-temperature calibration and averaging the local ratio 
T]/s over the trap size. In the figure we insert also the bound from string theory 
(dot-dot-dashed line), the low-temperature prediction of Rupak and Schafer 5 (dot- 
dashed line), and the high-temperature prediction of Bruun et al. 6 (dotted line). 
We plot also the results obtained with our model, Eq. d28b with Eq. (114b . for A = 
0.25 (solid line) and A = (dashed line). The figure shows that our model is in 
good qualitative agreement with the experimental data up to T/Tp c± 0.4. Both 
with A = 0.25 (solid line) and A = (dashed line) our model shows a minimum 
for r\/s ~ 0.44 at T/Tp ~ 0.27. Notice that the solid curve (A = 0.25) gives a 
reasonable agreement up to T/Tp ~ 0.9. 

We observe that the curve of X]/s vs T/Tp obtained by How and LeClair,- 
on the the basis of their version of Eq. ( 128b but with a very different procedure 
to calculate the scaled free energy <P(x), does not seem to increase as T/Tp goes 
to zero. Actually, a very recent calculation 31 of the shear viscosity from current- 
current correlation functions suggests that r\ /s at low T becomes small rather than 
exibiting the upturn. Nevertheless, the obtained theoretical values 31 of T\/s appear 
quite large with respect to the experimental ones. 



5 Conclusions 

We have described the elementary excitations of the unitary Fermi gas as made 
of collective bosonic excitations and fermionic single-particle ones. We have ob- 
tained an analytical expression for the Helmholtz free energy, showing that it is 
reliable to study the low-temperature thermodynamics of the unitary Fermi sys- 
tem up the critical temperature of the superfiuid phase transition. By using this 
free energy and simple scaling arguments we have derived the viscosity-entropy 
ratio r|/j as a function of the scaled temperature T/Tp. Contrary to other pre- 
dictions, our curve of r\/s vs T/Tp is in reasonable agreement with the available 
experimental data. 
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